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ABSTRACT 

Recently a new transport theory of cosmic rays in magnetized space plasmas 
extending the quasilinear approximation to the particle orbit has been developed 
for the case of an axisymmetric incompressible magnetic turbulence. Here we 
generalize the approach to the important physical case of a compressible plasma. 
As previously obtained in the case of an incompressible plasma we allow arbi- 
trary gyrophase deviations from the unperturbed spiral orbits in the uniform 
magnetic field. For the case of quasi-stationary and spatially homogeneous mag- 
netic turbulence we derive in the small Larmor radius approximation gyro-phase 
averaged cosmic ray Fokker-Planck coefficients. Upper limits for the perpendic- 
ular and pitch-angle Fokker-Planck coefficients and for the perpendicular and 
parallel spatial diffusion coefficients are presented. 

Subject headings: cosmic rays - diffusion - magnetic fields - plasmas - turbulence 
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1. Introduction 

The study of the cosmic ray transport in turbulent magnetic fields is crucial in many 
aspects of high energy astrophysics, such as the efficiency of cosmic ray diffusive shock 
acceleration, the modulation and penetration of low energy cosmic rays in the heliosphere 
and their confinement and escape from the Galaxy. 

A new theory of cosmic ray transport in magnetized plasmas extending the quasilinear 
approximation to the particle orbit has been recently published by one of us (Schlickeiser, 
2011) ( hereafter Paper 1 ). In Paper 1 the transport parameters of energetic charged 
particles in turbulent magnetized cosmic plasmas were derived for the case of an 
incompressible plasma, i.e. plasmas for which the component of the magnetic turbulence, 

— * 

SBz — 0, parallel to the guide magnetic field, Bq — BqCz, is set to zero. Here we present the 
generalization of the theory to the case of compressible magnetic turbulence with SB^ ^ 0. 

In Section 2 we briefiy review the theory developed in Paper 1. In Section 3 we 
obtain the gyro-phase averaged cosmic ray Fokker-Planck coefficients for a quasi-stationary, 
spatially homogenous turbulence under a Corrsin-type assumption on the nature of 
generalized orbits (Corrsin 1959, Salu and Montgomery, McComb 1990). Simplified 
formulas for the gyro-phase averaged cosmic ray Fokker-Planck coefficients are obtained 
in Section 4 assuming that the magnetic turbulence is asymmetric, while the quasilinear 
limit of the coefficients is shown in the Appendix. In Section 5 from the Fokker-Planck 
coefficients we derive upper and lower limits for the perpendicular and parallel spatial 
diffusion coefficients. In Section 6 we compare the relative importance of mirror forces and 
turbulent scattering for the cosmic ray transport in interstellar plasmas. 
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2. The gyro-averaged Fokker-Planck equations 

2.1. Equations of motion of a particle in magnetic fields 

For the following treatment we remind shortly the equation of motion of charged 
particles of mass m, charge g, and Lorentz factor 7 = a/1 + {p/mcY in a uniform guide 
magnetic field Bq = Bqc^ = {0,0, Bq). A random magnetic field, 6B, is superposed to the 
guide field. 
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The scalar product of Eq. ([T]) with p readily yields p = \p\ =const., v =const. and 
7 =const.. Introducing the constant relativistic gyrofrequency = and scaling the 
turbulent fields in units oi Bq, 5h = -§7 we obtain 
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For the time evolution of the particle pitch-angle cosine /i 
a.Tcta.n{vy/vx) this implies 
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and 



^ = — + h^{t), h^{t) = —Vtdbz H ^ (cos 6bx + sin 6by) . (4) 

dt ^1 — ^2 

with the two random forces hfj_{t) and h^{t). 
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In the coordinates of the guiding center 



-Vx 

\ / 



Eqs. ([2]) become 
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dt 
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Indicating = [X, Y] with i,j = 1, 2, Eq.(l6]) provide the two additional random force 
terms hi{t), proportional to the turbulent magnetic field components 



dXi 



hi{t) = v,{t)6bi{t) - v,{t)6b,{t), 



(7) 



2.2. The ensemble-averaged particle distribution function 

The description of the cosmic ray transport within a large-scale guide magnetic field, 
which is uniform on the scales of the cosmic ray particles gyradii Rl = v/\Q\, is given by 
the solution of the Vlasov (collision-free Boltzmann) equation for the particle distribution 
function F (Hall and Sturrock 1968, Schlickeiser 2002). In spherical momentum coordinates 
{X,Y, z,p, fi,(f)) the Vlasov equation reads (Hall and Sturrock 1968, Achatz et al. 1991) 



OF OF OF d d 

'dt^ ^/i^-f^^+p"'^[ApWi^]+^[/^aWi^]-Qo(^,X,F,p,/i,0,t) = O, (8) 
where i/a € [X, Y, fi,(j)] and 
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Qo{z, X, F, p, /i, 0, t) = 5o(2;, X, F, p, /x, 0, t) - Af^F - TZqF 



(9) 



accounts for sources and sinks (5*0) and the effects of the mirror force (Aq) and momentum 
loss processes (T^o); where the latter two operate on much longer spatial and time scales 
than the particle interactions with the stochastic fields. In Eq. ([8]) we use the Einstein sum 
convention for indices and the short notation = (d/dx^). x^^a £ [fi',P,X,y] represent the 
four phase space variables fi,p,X,Y with non- vanishing stochastic fields h^{t). 

The particle distribution function, F{X,Y, z,p, fi,(f)), varies in a irregular way under 
the influence of the stochastically fluctuating fields, hy{t). However, we do not look for 
the detailed function F, but rather we look for an ensemble-averaged solution, < F >, an 
expectation value of of Equ. [H which results from averaging over different realizations of 
the fields hy{t) with the same statistical averages. In the following treatment we will keep 
only first-order terms in the fluctuating quantities 5F 



and in the turbulent fields, hi^{t). In other words we will consider the case of weak 
turbulence or quasilinear approximation. Moreover we will also neglect electric fields 
{hp = 0). As shown in details in Paper 1, under the assumption of weak turbulence the 
ensemble- averaged solution, < F >, can be obtained by solving the kinetic equation 



F=<F> +SF 



(10) 



h<F> +v^id, <F> -Qd^ <F> -Qq{z, X, Y,p, /i, 0, t) 
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(11) 



where Xj, G [X, Y, /i] and the Fokker-Planck coefficients are given as 




(12) 
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The time-integration operator in Equ. [12] is performed over a generalization of the 
unperturbed gyrocenter orbit in the uniform magnetic field with deviations of the gyrophase 
given by 

Xs = X, Ys = Y, Zs = Z + v^{s -t), Ps = p, /is = /i, 0s = - ^{s - t) + 6(/)(t - s), (13) 

that contains the additional arbitrary gyrophase variation 6(f){t — s), with 6(f) = for s = t. 
Fourier transforming the stochastic force in space, the time integral in Equ. [T2] becomes 



t f rt 

3; 



dshfj{s)= I d k / duH„{k,s) 

to J J to 



X exp 



ik ■ X + iv^k\\ {s — t) + ik±v \/l — jj,"^ / dw cos — (f) + ri{w — t) — Sep {t — w)) 



(14) 



where the particle position is given as 



xis) 



\ 



fX + v^l-fi'^ fdw cos{cf> ~ n{w - t) + 6(f){t - w))\ 
Y + v^yl- fi^ /' dw sin(0 - n{w - t) + 5(p{t - w)) 
Z + vfi{s — t) 

and where we have introduced cylindrical coordinates for the wavenumber vector 
k = {k± cosip, k± sinip, k\\) and the particle velocity. 



(15) 



As explained in details in Paper 1, in the small Larmor radius approximation (Chew et 
al. 1956, Kennel and Engelmann 1962) the distribution functions are independent of to 
lowest order and can then be expanded as 



(16) 



and the Larmor-phase-averaged equation becomes 
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(9 (9P 

dtFo + vfid^Fo - Q(z, X, F, p, /i, t) = --^D^,-^ - (17) 



with the gyro-averaged source term 



1 Z"^'' 

Q{z,X,Y,p,fi,t) = — d<PQo{z,X,Y,p,fi,<p,t), (18) 
and the gyro-averaged Fokker-Planck coefficients 



Daa = 3f?7^ / d(t)Paa = ^7^ d(j)< K{t) / ds K{s) >, (19) 

Jo Jo Jto 

where we replaced h^{t) = h*^{t) by its complex conjugate because the stochastic forces 
are real-valued quantities. 

The generalization of the time integral in Eq. [12] from the unperturbed motion of 
the gyrocenters in the guide magnetic field to arbitrary gyrophase motions of particles is 
possible essentially because of the gyrophase-averaging in Eq. [T3 As demonstrated in Paper 
1 the considered general particle gyrophase motion then only modifies the arguments of 
trigonometric and Bessel functions as compared to the quasilinear approximation of particle 
orbits. 



3. Derivation of Fokker-Planck coefficients for compressible magnetic 

turbulence 

Following the approach of Paper 1 we now derive the Fokker-Planck coefficients for the 
case of compressible magnetic turbulence. We make the following assumptions on the nature 
of the turbulence: the turbulence is quasi-stationary, meaning that the correlation function 
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< hl(t)ha{s) > depends only on the absolute value of the time difference |t — s| = |r|, so 
that with the substitution s = t — r we find for Eq. ([T^ 



= SR— / d(P ds < K{t)hl{s) >= 3fJ— d(P ds < K{t)hl{t-T) > (20) 

27r Jo Jto 27r Jq Jq 

As second assumption we use that the turbulent magnetic fields are homogenously 
distributed, meaning that independent from the actual position of the gyrocenter at time 
t the particles are subject to turbulence realizations with equal statistical properties. This 
allows us to average the Fokker-Planck coefficients over the spatial position of the guiding 
center using 

d^X e'^^-^^-^ = 5{^! -k), (21) 



(27r)3 

implying that turbulence fields at different wavevectors are uncorrelated. As explained in 
details in Paper 1, the Fokker-Planck coefficients then become 



X exp 



D^^ = ^— [ d(f) [ d^k [ ' dT < H^(k,t)Hl{k,t - T)e'''^'''\\^ 
2vr Jo J Jo 

-ik^v^l- li^ (j dw cos(^ - + Q{w - t) - S(f){t - w)) + ^^^^^^ 



> 

(22) 



A third assumption concerns the nature of the particle orbits, i.e. we consider only 
orbits where 6(j){w) does not depend upon the fluctuating fields, so that the ensemble 
averaging in Eq. fl22|) involves only the 2nd order correlation functions of the stochastic 
fields. This is generally called the Corrsin independence hypothesis (Corrsin 1959, Salu and 
Montgomery, McComb 1990). 

With ^ = t — w and the abbreviation 

G{0 = n^ + Sm (23) 
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the Fokker-Planck coefficients fl22l) then are 



'2tt 



2tt 



t-to 



dr < H^{k,t)H;{k,t - t) > 



X exp 



tvjj^kiiT + tkj_v\/l — fi'^ I / d^ cos{(j) — ip + G{^)) 



sin(0 — ip) 



(24) 



Later we will also assume that the turbulence has a finite decorrelation time such 
that the correlation functions < h^{t)h'^(t — t) >— t- fall to a negligible magnitude for 
r — )■ cxD, so that the upper integration boundary in the r-integral can be replaced by infinity 



D,. 



/*2'jr poo 

— d(j) dr < K{t)hl{t - r) > . 
Jq Jo 



(25) 



'0 ^0 

We remark that diffusive transport of cosmic rays happens if the turbulence is 
quasi- stationary and has a finite decorrelation time tc, because the resulting gyro-averaged 
Fokker-Planck coefficients at large times t — ^ t^. no longer depend on the time-difference 
t-to. 

The equations of motion of the guiding center (Eqs. [7j) can be written as 



dXj 
~df 



hi{t) = vfi5bi{t) - Vi{t)5bz{t), 



if 



(26) 



Vi{t) = v\/\- y? cos{{i - 1)- - 0) 
where i = [1,2], whereas the pitch- angle random force (Eq. ([3])) is 



(27) 



^ = h^{t) = n^/T^ (cos (P5b2it) - sin (l)5b^{t)) , 



{2i 
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3.1. Individual gyro-averaged Fokker-Planck coefficients 

The Fourier transforms of tlie stocliastic fields in (l26l) and fl28l) are 

Hi{k,t) = vfibi{k,t) - Vi{t)bz{k,t), 

H*{k,t-T) = vfj.b*{k,t - r) - Vi{t - r)bl{k,t - t), 
Hf,{k,t) = fiv^l -/i2 (cos{(f))b2ik,t) - sm{(j))bi{k,t) 



H;{k,t-r) 
where 



- fi^ cos(0 + G{T))b;{k, t-T)~ sm{(f) + G{T))bl{k, t-T)) (29) 



(t) = v^l-fx^ cos - 1)| - 0) , Viit - r) = t;v/l-/i2 cos ((z - 1)^ - (0 + G(r))) 



In terms of the magnetic field correlation tensor 



<bi{k,t)b*{k,t-T) >=P,j{k,T) 

we then obtain for the perpendicular Fokker-Planck coefficients 



(30) 



= SR— / d^k 



t-to 



dr e 



IVflku T 



2tt 



d(j) Hij{k,T) 



X exp 



zk^v^/l - /i2 ( / di cos(0 - ^ + G(0) 



sin(0 — ?/;) 



where 



(31) 



H,,{k, r) = ^I'P^Jik, r) - Vl - cos ( (j - 1)- - (0 + G(r)) ) Pi,(A;, r) + 



7r 



-/X Vl - /i^COS ( (t - 1)- - 0) P,jik,T) + 



TT 



+ (1 - /i^) COS (z - 1)- - COS (j - 1)- - (0 + G(r)) P,,(A;, r) (32) 
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The mixed Fokker-Planck coefficients are instead given as 



where 



and 



where 



D. 



2-K 



X exp 



Hi^{k,T) = iicos{(p + G{t)) Pi2{k,T) - II sin {(P + G{t)) Pii{k,T) + 
- yr^cos ((i - 1)| - 0) cos((/) + Gir))P,2ik, r) + 
+ Vl-/^'cos ((i - 1)| - 0) sm{(f> + G{r))P,i{k, r) 



27r 



f27r 



X exp 



sin(0 — ■?/^) 



— * — * — * 

H^ii{k, t) ^ n cos (t>P2i{k, t) - n sin (f)Pu{k, r) + 
- V^T^COS 0COS ((i - 1)| - (0 + G'(t))) P2.(^, t) + 
+ yr^sin^cos ((^ - 1)^ - (0 + G'(r))) Pu(^,r) 



(33) 



(34) 



(35) 



The Fokker-Planck coefficients parallel to the direction of the guide magnetic field are 



/rt—to r2n 
d'k / dre'^'^'^w^ / d<f>H^^{k,T) 
Jo Jo 



X exp 



ikxv^/l^ (J d^ cos(0 - V + G'(O) - 



sin(0 — ip) 



(36) 



where 



H^^,{k, t) = cos cos(0 + G{T))P22{k, t) + sin sin(0 + G'(t))Pii(^, t) 



- sin (f) cos(0 + G{T))Pi2{k, r) - cos sin(0 + G{T))P2i{k, r)) 

The 0-integrals are calculated in Appendix A of Paper 1 and yield 

,2 r rt-to 



dre'''^''\\^Iij{k,r) 



(37) 







where 



Iij{k, t) = [2fj,^Pijik, t) + (-1)^1^-^1(1 - /x^) cos (l i - J I I - G'(t)) P,,(^, r)] Jo(Z) + 
-2zVl-//2//[(-l)'~'sin (^(^ - 1)| - (V' + arcsin(^))^ P.,-(^,t)+ 



(-l)^'-'sin I (j - 1)^ - (V. + G(t) +arcsin(^)) ) P,,(^, r)] Ji(Z) + 



TT 



(2^^ + G{t) + arcsin(^)) ) J2(Z) 



in 



(38) 
(39) 



where 



Iilx{k,T) = Vl - 



sin ((^ - 1)| - G(r)) P,i(^,t) + (-l)^cos ({i - 1)| - G(t)) P,2(^,t)] Jo(Z) 

Ji(^) + 



sin ( + G{t) + arcsin(^) j Pi2(^, r) + cos + G{t) + arcsin(^) ) Pii(/i:, r 



(-1)'-^ sin ( (i - 1)| - (2^^ + G{t) + arcsin(^)) ) P,,{k, r) + 



- cos ( (i - 1)| - (2^^ + G{t) + arcsin(^)) ) P,2(A;, r) | J2(Z) 



(40) 



v|J,fl^/l — ^j? 



t-to 

d'k I dre'''^'''\\^Ii,i(k,T) 





(41) 
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where 



and 



-l)^sin ((i - 1)| - G(t)) Pi.(^,t) - cos ((i - 1)| - G'(t)) P2.fiT)] Jo(Z) 

]Ji(Z) + 

v/r^[(-l)^-^ sin (^{i - 1)| - (2^^ + G(t) + arcsin(^))^ P^k, r) + 



sin ( + arcsin(^) j P2i{k, r) + cos (ip + arcsin(^) ) Pii{k, r) 



- cos ( (i - 1)| - {2^ + G{t) + arcsin(^)) ) P^^ik, r) \ J2{Z) 



X 



cos{G{t))Jo{Z) - cos(2V' + G{t) + 2 arcsin(^)) J2(Z)^ Pn(^, r) 

+ ( cos(G'(t)) Jo(Z) + cos(2V' + G{t) + 2 arcsin(^)) J2(Z)^ P22(^, r) 

sin(G'(r)) Jo(Z) + sin(2V' + G{t) + 2 arcsin(^)) J2(Z)^ P2i(^, r) 

+ (^sin(G'(T)) Jo(Z) - sin(2V' + G{t) + 2 arcsin(^)) J2(Z)^ Pi2(^, r) 
respectively, where JniZ) denotes the Bessel function of the first kind and order n, 



(42) 



(43) 



Zi = k±v^/l - / cos (G'(O) 



and 



(44) 



Z = k±v\/l — 



di cos(G(0) 



1 



di sin(G(0) 



1/2 



(45) 
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4. Axisymmetric turbulence 



Useful formulas can be obtained by assuming that the turbulence is asymmetric, 
meaning Pa^(/c, r) are independent of the wave phase ip 



Pap{k,T) = P«^(A;||,/i;_L,r), 



(46) 



The integration over ip of the general formulas ( |38l) - ( l43l) then provides 

rt-to 



OO POO 



dk\\ I dk± /c. 

'o ^0 



X 



TT 



2/x^P.,(fc||, r) + (-1^1-^1(1 - ;,^) cos ^ - j | - - G{r) P..{k\\M. r) , (47) 



OO />oo 



-OO Jo 



dre^^^'^ii" Jo(Z) 



sin ( (^ - 1)| - G(r)) P.i(A;|i, fc^, r) + (-1)^ cos ((z - 1)| - G(r) ) P.2{k\\,k^, r) (48) 



OO /"OO 



D., = ^7rvn{l- 12^) dk\\ dk^k^ dr e''"'^^\^ Jo{Z) 



t-to 



-OO JO JO 



-1)' sin ( (z - 1)- - G'(r) j Pi,(A;i,, A;^, r) - cos - l)^ - G{r)j P2,{k\\,k^, r) (49) 



and 



OO poo 



t-to 



D^^ = ^7cv'Q'\l-fi'') I dk\\ I dk±k^ I tire^^^^ii" Jo(Z) 

J~oo Jo Jo 

cos(G(r)) (Pii(A;,|, A;^, r) + P22(A;||, A;^, r)) + sin(G(r)) (Pi2(A;||, A;^, r) - P2i(A;||, A;^, r)) 

(50) 



Introducing the left-handed and right-handed polarized stochastic magnetic field 
components 
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5hL,R = [6bi ± i5b2] 



so that 



(51) 



2Pll = Pu + P22 + iP2i - iPi2, 2Prr = Pu + P22 + iPi2 - ^P. 



21, 



(52) 



we obtain for the pitch-angle Fokker-Planck coefficient ( 150 



00 poo 



00 JO 



t-to 



X 



g,(.Mir+GWp^^(^l^^ r) + e<'"^^ii^-^MpRR(A;||, fcx, r) 



(53) 



5. Upper and lower limits of the general Fokker-Planck coefficients in the 

diffusion limit 

If we now consider a magnetic field fluctuation decorrelation time tc = (Schlickeiser 
and Achatz 1993, Bieber et al. 1994) 

P.,{k,T) = P^^{k)e-^\ (54) 

then in the diffusion limit t — to ^ tc the general Fokker-Planck coefficients ( H7j) - ( |53l) in 
asymmetric turbulence become 



00 POO 



D^^ = nQ\l - fi') dk\\ / dk^k^ dTJo{Z)e-"'^ 



X 



cos{v^k\\T + G{T))Plj^{k\\, k±) + cos{vfik\\r — G{T))P^j^{k\\, fcj 



(55) 



and 
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oo Jo Jo 



X 



1 



(56) 



/OO /'OO /"OO 

dA;|| / dk^k^ I rfr Jo(^)e™^'=ii""^" 
■oo Jo Jo 



(57) 



/OO /"OO /"OO 

oo Jo Jo 



X 



Note that in Eqs ( I55ll58l) we consider only the real part of the integral. 



(5^ 



Because of the existence of the finite turbulence decorrelation time 7"^, the correlation 
functions Pij{k^\, k±, r) fall to a negligible magnitude for r — )■ 00, allowing us to replace the 
upper integration boundary in the r-integrals in Eqs. f l55|) . fl56|) . f l57|) and f l58|) by infinity. 
We recover the diffusion limit which is valid for times t — to ^ 7^^. 

If Jo{A) < 1 and cos(x) < 1 Eqs. ([55]- [581) become 



D < D"'^'' = ^ ^ 



/oo /»00 /"OO 

dk\i / dk±k±P^j{k\\,k±) / (ir Jo(2') cos(t'/iA;|jr) 
-00 Jo Jo 



-'yr 



7 J-oo Jo 



+ 



27r / dki I dk^k±P^^{k\\,k±) 



dr Jo{Z) 



COS 



(f/iA;||r + G(t)) + cos {yjj,k\\T — G{t)) 
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and 



7 



+ 



27 



(59) 



dk±k± / dr Jq{Z) 



-7r 



[cos [v^k\\T + G{t)) + COS [v^k\\T - G{t))) 5i2-Pfi(A;||, /c±) 
\cos {yiJik\\T + G'(t)) + cos (f/i/cyr — G'(t))] (5jiP°2(^|h ^-l) 



27 



(60) 



^■Qfl _ m2\ /-oo /.oo /-cx) -7r 

D^, < = \ ^ ^ 27r / dk\\ / dfc^ fc^ / dr MZ) -— 

^ J~oo Jo Jo ^ 

[cos {yiik\\r + G'(r)) + cos {yiik\\r - G'(r))] 5i2-Pii(^lh ^-l) 
[cos (f/i/cyr + G'('r)) + cos [y^k\^T — G'('r))] (5iiP2z(^||) ^J- 
= 1^ \0i2b\.^ - kxbb^^\ 



(61) 



27 

^2(l-;,2^ 



/OO f'OO 
dk\\ / dfcx fcx [PlL{k\\,k,_) + P^^(fc||, A;j 
■00 -/o 



27 



where 



(62) 



//»00 /"OO 
t/3A;P;,(A;||, A;^) = 27r / rfA;,, / rfA;^ k^Pl,{k\\,k, 
J -00 Jo 



(63) 



According to the diffusion approximation (Schhckeiser 2002), neglecting the inffuence of 
the mirror force contribution Ao in Eq. (j8]), the perpendicular spatial diffusion coefficients 
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for the isotropic part of the cosmic ray phase space density are given by the pitch-angle 
average 



1 



dfiD,,{fi) (64) 
^ J-i 



From Eq. EH] we find the upper hmits 



K,, < K--^ = - {Sbl + SblA,) (65) 

The parallel spatial diffusion coefficient for the isotropic part of the cosmic ray phase 
space density is given by the pitch-angle average 

and its lower limit does not change with respect to the incompressible case examined 
in Paper 1. 



Kii > ^r-^ = — ^ (67) 



6. Mirror forces and turbulent scattering in the solar wind plasma 

Mirror forces are produced by large-scale spatial variations of the guide magnetic field. 
The perpendicular component of the mirror force generates gradients and curvature drifts 
of the cosmic ray guiding center ( Boyd and Sanderson 1969 ). In the presence of mirror 

(T) 

forces the diffusion coefficients are given by the sum of the turbulent contribution, k^iu , and 
of the contribution due to the mirror forces, kl!^\ 
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= + '^if ^ (68) 

We will now compare the effect of mirror forces and of the turbulent contribution, kj^^^^^, 
calculated in Sect. 5, on the properties of cosmic ray transport in the solar wind plasma. 

For mirror forces Schlickeiser and Jenko (2010) showed that in the case of a symmetric 
choice of the pitch-angle Fokker-Planck coefficients the ratio of the perpendicular mirror 
spatial diffusion coefficient to the parallel turbulent spatial diffusion coefficient is given by 
the derivatives of the cosmic-ray Larmor radius. In particular, considering the case of a 
magnetic power spectrum of Alfvenic slab turbulence P{k) oc with s < 2, Schlickeiser 
and Jenko 2010 obtained that the ratios of the non-zero perpendicular to parallel spatial 
diffusion coefficients are 



(T) Q — s SLo 



) 



^YY 2 S ^ i?£ ^ 



'^zz 



(69) 



where we introduce the perpendicular magnetic field scale lengths (Schlickeiser and 
Jenko 2010) 

ox 

Lr' = -B-'f (70) 
Summing the diagonal terms of the diffusion matrix we obtain 



AM) , AM) 2 - s /?r 2 / 1 1 



6-.^ 3 ' \Ll LI 
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where we remind that the cosmic ray gyroradius is defined as 



From Eq. [7T] we have 

3 ' 6-s V^f L 



smce 



The relevant magnetic field random irregularities for the cosmic ray transport properties 
are the fast magnetosonic waves ( Lee and Volk 1975, Cho and Lazarian, 2003 ). If we 
consider isotropic magnetosonic waves ( Schhckeiser 2002 ) 



Prr = ,l cOs9 



P =^ 
87rA;2 

P^^ = li§sine2 (75) 



then the random contributions to the field irregularities are 



^bl^ = \j dfxfi"^ j dkg{k) = ^ j dkg{k) 

^Ky = \j j dkg{k) = ^ j dkg{k) 
Sbl = \ dii (1 - /i^) j dkg{k) = dkg{k) 



(76) 



Using the upper and lower limits for the diffusion coefficients from random turbulent 
forces in Eqs. ( 165|) and (!67|) and applying it for the case of fast magnetosonic waves we have 



- 22 - 



(J,) mm / (T)max (-p)max 



> 



,vR, 



1 + 



256' 



,2 

22 



" XX ' yy 



(77) 



Taking the ratio of Eq. [73] with Eq. [77] we obtain a relation for the product of 
perpendicular diffusion coefficients independent of kzz 



+ K 



(M) 
YY 



max 

XX 



> 



3 ' 
2 - s 



2(6 



18(6 — s) min[Li,L2 



(7J 



Perpendicular spatial diffusion is thus dominated by turbulent forces at low particle 
momenta, where the gyroradius is less than the minimum of the perpendicular magnetic 
field focusing lengths. Alternatively, at high momenta, where the gyroradius is larger than 
the minimum of the perpendicular magnetic fierld focusing lengths, perpendicular diffusion 
is dominated by the mirror force contribution. 



7. Summary and conclusions 

In a large-scale magnetized plasma the description of cosmic ray transport is given by 
the solution of the Vlasov equation for the particle distribution function. The influence of 
stocastically fluctuating fields on the particle distribution function can be studied by looking 
for an ensemble-averaged solution of the Vlasov equation, which results from averaging over 
different realizations of turbulent fields with the same statistical properties. 

In the small Larmor approximation it was shown in Paper 1 that one can obtain the 
solution of the Vlasov equation for arbitrary gyrophase motions of the particles, extending 
the quasilinear approximation to the particle orbit. In Paper 1 the transport parameters 
of energetic charged particles in turbulent magnetized cosmic plasmas were derived for the 
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case of an incompressible plasma, i.e. plasmas for which the component of the magnetic 
turbulence, SB^ = 0, parallel to the guide magnetic field, Bq = Bqc^, is set to zero. 
Here we present the generalization of the theory to the case of compressible magnetic 
turbulence with 6Bz 7^ 0. Under the assumption that the turbulence is quasi-stationary and 
homogeneous we have obtained the gyro-averaged Fokker-Planck coefficients for a Corrsin 
type of generalized orbits. For an axisymmetric turbulence we have derived upper and lower 
limits for the perpendicular and pitch-angle Fokker-Planck coefficients in the diffusion limit. 
We have shown upper and lower limits for the perpendicular and parallel spatial diffusion 
coefficients, respectively, describing the spatial diffusion of the isotropic part of the cosmic 
ray phase space density. Finally using the upper and lower limits for the turbulent motion 
we compare the effects on the transport of cosmic ray particles of the turbulent and of 
mirror forces, if the latter cannot be neglected. 

This work was partially supported by the German Ministry for Education and Research 
(BMBF) through Verbundforschung Astroteilchenphysik grant 05A11PC1 and the Deutsche 
Forschungsgemeinschaft through grant Schl 201/23-1. 



8. Appendix: Quasilinear limit 



Following the approach in Paper 1 and assuming 50 = for the particle orbit in Eq. 



G(0 = n^ + sm 



(79) 



we obtain the quasilinear approximation to the particle orbits (Shalchi and Schlickeiser 
2004). The argument of the Bessel functions of first kind 



Z = k±v\/l — fi"^ 



cos(G(e)) 



1 



sin(G(0) 



1/2 
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at order n = becomes 



Z = Zn 



sin (fir) + (1 -cos2(fir))- 



nl/2 



[2 (1 - cos(fir))] 



1/2 



U """^T^ 



The perpendicular Fokker-Planck coefficients (H7|) then become 
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The mixed Fokker-Planck coefficients (HHlandHOl) are 
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sin ( (z - 1)| - fir) P,i(A;j|, A;^, r) + (-1)' cos {{% - 1)^ - ^r) P,2(A;||, fc^, r) (83) 
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-1)' sin ( (z - 1)| - fir) Pi.(A;|j, fcx, r) - cos ((z - 1)| - fir) P2.(A;||, A;^, r) 



M) 



whereas the Fokker-Planck coefficient (l53l) reduces to 



/■t-to /-oo /-oo /9i',7i, /Or 
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Note that the quasilinear Fokker-Planck coefficients in axisymmetric turbulence no 
longer involve infinite sums of products of Bessel functions which enormously facilitates 
their numerical computation for specified turbulence field correlation tensors. 



We now explicitly calculate the different Fokker-Planck coefficients 
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Using the Bessel function addition theorem (see Appendix in Paper 1) with ri 

k±v±/fl and 9 = fir then 



j.(z„) = j„ (^ [2(1 - cosjix)]"^) = f: (^) 
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SO that the perpendicular Fokker-Planck coefficients can be written as 
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n=— oo ^ 

E e^'^ii^ir-K _ j.^^ ^^^^ P..(fc||,fc„r)], (97) 
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Using the addition theorem flMj) the pitch angle Fokker-Planck coefficient f l85|) becomes 



°° /"OO /"CXD /-C 

Dgf = 3?7^^]2(l_^2^ ^ / c/A;|| / dk^k^ 

„-_^J-oo Jo Jo 



(99) 



The pitch angle coefficient in fl99|) does not change with respect to the case of incompressible 
plasma treated in Paper 1. As remarked in Paper 1 the pitch angle coefficient agrees exactly 
with Eq. (12.2.5) of Schlickeiser (2002) for incompressible, axisymmetric turbulence. 



-28- 

REFERENCES 

Achatz, U., Steinacker, J., Schlickeiser, R., 1991, A& A 250, 266 

Bieber, J. W., Matthaeus, W. H., Smith, C. W., Wanner, W., Kallenrode, M.-B., 
Wibberenz, G., 1994, ApJ 420, 294 

Burger, R. A., Hattingh, M., 1998, ApJ 505, 244 

Chew, G. F., Goldberger, M. L., Low, F. E., 1956, Proc. Roy. Soc. (London) A236, 112 
Cho, Jungyeon; Lazarian, A., 2003, MNRAS 345, 325C 

Corrsin, S., 1959, in: Atmospheric Diffusion and Air Polution, Advances in Geophys., Vol. 
6, eds. F. Frenkiel and P. Shepperd, Academic Press, New York, p. 161 

Hall D. E. and Sturrock, P. A., 1968, Phys. Fluids 10, 2620 

Jaekel, U. and Schlickeiser, R., 1992, J. Phys. G:Nucl. part. Phys. 18, 1089 

Kennel, C. F., Engelmann, F., 1966, Phys. Fluids 9, 2377 

Lee, M. A. and Volk, H. J., 1975, ApJ 198, 485 

McComb, W. D., 1990, The Physics of Fluid Turbulence, Oxford Univ. Press, New York 

Salu, Y., Montgomery, D. C., 1977, Phys. Fluids 20, 1 

Schlickeiser, R., 2002, Cosmic Ray Astrophysics, Springer, Berlin 

Schlickeiser, R., Achatz, U., 1993, J. Plasma Phys. 49, 63 

Schlickeiser, R., Jenko, F., 2010, J. Plasma Phys. 76, 317 

Schlickeiser, R., 2011, ApJ 732, 96 



-29- 

Shalchi, A., Schlickeiser, R., 2004, A& A 420, 821 



This manuscript was prepared with the A AS macros v5.2. 



